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ABSTRACT 

We constrain Galactic foreground contamination of the Python V cosmic microwave 
background anisotropy data by cross correlating it with foreground contaminant emis- 
sion templates. To model foreground emission we use 100 and 12 /im dust templates 
and two point source templates based on the PMN survey. The analysis takes account 
of inter-modulation correlations in 8 modulations of the data that are sensitive to a 
large range of angular scales and also densely sample a large area of sky. As a conse- 
quence the analysis here is highly constraining. We find little evidence for foreground 
contamination in a analysis of the whole data set. However, there is indication that 
foregrounds are present in the data from the larger-angular-scale modulations of those 
Python V fields that overlap the region scanned earlier by the UCSB South Pole 1994 
experiment. This is an independent consistency cross-check of findings from the South 
Pole 1994 data. 

Subject headings: cosmology: observation — cosmic microwave background — diffuse 
radiation — dust, extinction 



1. Introduction 

While cosmic microwave background (CMB) anisotropy data have started to provide interesting 
constraints on cosmological parameters (see, e.g., Podariu et al. 2001; Page 2002; Mukherjee et 
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al. 2002b, 2003; Benoit et al. 2003; Ruhl et al. 2002; Kuo et al. 2002; Slosar et al. 2002 for 
recent results, and Peebles &: Ratra 2003 for a review), Galactic emission foreground contaminants 
in them are still not well understood. Robust constraints on cosmological parameters from these 
data require a better understanding of the effect of these contaminants. 

In this paper we study foreground contaminants in the Python V (hereafter PyV) CMB 
anisotropy data. Python V is the latest of the Python experiments at the South Pole. Coble 
et al. (1999, 2003) describe the PyV experiment, observations, and data reduction. Dragovan et al. 
(1994), Ruhl et al. (1995), and Piatt et al. (1997) describe Python I III and Rocha et al. (1999) 
and Mukherjee et al. (2003) derive constraints on cosmological parameters from these data. Coble 
et al. (2003) also describe the procedure used to create maps of the sky with PyV and Python III 
data; these maps were compared to infer consistency and indirectly deduce the lack of significant 
foreground contamination in these data. 

The PyV data are acquired at a frequency of ~ 40 GHz. Two regions of sky covering 598 deg 2 
in the southern hemisphere were observed (the "main PyV" region and a smaller region, the fields 
labelled 'sa', 'sb', and 'sc' in Coble et al. 2003, that encompasses the region scanned earlier by 
the UCSB South Pole 1994 experiment — hereafter the "SP94 overlap" region). 690 fields were 
scanned in all (345 effective fields were scanned with 2 detector feeds separated by 2?8 in azimuth 
on the sky) with an asymmetric Gaussian beam of FWHM 0?94 x 1?02. Once the telescope was 
positioned on each field, the chopper smoothly scanned in azimuth with a throw of 17°, 128 samples 
were recorded in each chopper cycle, and 164 chopper cycles of data were taken of a given field. 
The densely sampled data were then modulated in software using the first eight cosine harmonics of 
the chopper cycle (hereafter modulations — or in Tables mods — 1-8). The modulations approach 
used has the advantage of filtering out some of the contaminants in the time stream, and also 
provides a rapid means of compressing a large amount of data into a more manageable size. All the 
modulations, other than the first, were apodized by a Hann window to reduce ringing in multipole 
space and down weight data taken during chopper turnaround. The resulting data are sensitive to 
angular scales corresponding to multipole moments ranging from I ~ 40 to I ~ 260. More details 
about the particular observing strategy employed and results found, such as the angular power 
spectrum of the data, may be found in Coble et al. (2003); see Souradeep & Ratra (2001) for 
details about the window functions. 

For this foregrounds analysis, we follow the general method outlined for example in Hamilton 
& Ganga (2001) and Mukherjee et al. (2002a), doing a multi-modulation analysis here rather 
than the multi- (frequency) channel analysis discussed in these papers. In doing so we extend the 
preliminary estimate presented in Coble et al. (1999). Here we use the technique of Souradeep 
& Ratra (2001) to account for correlations between modulations, in data that are sensitive to a 
substantial range of angular scales, as well as compare to what has been found from the SP94 
experiment (Gundersen et al. 1995; Ganga et al. 1997) data about foreground contamination in 
part of the PyV region (Hamilton & Ganga 2001; Mukherjee et al. 2002a). 
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CMB data have previously been correlated with foreground templates (DMR: Kogut et al 
1996a, 1996b; 19GHz: de Oliveira-Costa et al. 1998; Saskatoon: de Oliveira-Costa et al. 1997; 
OVRO: Leitch et al. 1997, Mukherjee et al. 2002a; SP94: Hamilton & Ganga 2001, Mukherjee 
et al. 2002a; Tenerife: de Oliveira-Costa et al. 1999, 2002, Mukherjee et al. 2001; QMAP: de 
Oliveira-Costa et al. 2000; MAX: Lim et al. 1996, Ganga et al. 1998; and Boomerang: Masi et 
al. 2001) 8 . Correlations between CMB data and infra-red emission seem to be roughly consistent 
with free-free emission, spectrally, over a wide range of frequencies (10 to 90 GHz) and angular 
scales (7° to 7'), with some evidence for a contribution from spinning dust emission with a peak 
around 15 to 20 GHz. In general though the contamination from foregrounds has not been found 
to be large in any experiment, but on the detail level residual foregrounds can cause problems with 
parameter estimation and non-Gaussianity tests for high precision CMB data. 9 

To model diffuse Galactic emission, we use the 100 /im IRAS+DIRBE map (Schlegel, Finkbei- 
ner, & Davis 1998) as a tracer of thermal emission from interstellar dust, and the 12 ^m map 
(D. Finkbeiner, private communication, 2000) as a tracer of emission from ultra-small dust grains. 
Emission from small dust grains is still under study (Finkbeiner et al. 2002). Such grains may con- 
tribute significantly at microwave frequencies according to a model by Draine & Lazarian ( 1998a, b), 
and this has been reviewed from the CMB data point of view by Kogut (1999) and Draine & Lazar- 
ian (1999). The Haslam 408 MHz map (Haslam et al. 1981) of synchrotron emission was not used 
as it does not have enough resolution for all the modulations of the PyV experiment to be simulated 
on it. At the same time synchrotron emission is not likely to contribute significantly at 40 GHz. 
We have tried to use the SHASSA H a data (Gaustad et al. 2001) as a tracer of free-free emission, 
but found that the template has insufficient resolution for our purpose. 

We also use two point source templates created from the PMN survey (Wright et al. 1994). 
The one called PMN has been converted to 5Tqmb (equivalent temperature fluctuations in the 
CMB at 40 GHz) using the spectral indices given in the survey, while the one called PMN0 is 
converted to a flux at 40 GHz assuming a flat spectrum extrapolated from the flux measurement 
at 4.85 GHz. The assumption of a flat spectrum is conservative in that it is likely to overestimate 
the flux at 40 GHz. Neither case is correct, since spectral indices have not been measured for all 
of the sources, in which case a flat spectrum is assumed. The two cases cover a reasonable range 



8 DASI (Halverson et al. 2002) and VSA (Taylor et al. 2002) are interferometric experiments, that observed at 
frequencies of 26-36 GHz and are sensitive to multipoles of ~ 100 to 900. Point sources are the dominant source of 
contamination and the effect of contamination from diffuse Galactic foregrounds upon the CMB power spectrum is 
inferred to be small for both data sets. A full cross-correlation analysis of VSA data is underway (Dickinson et al., 
in preparation). 

The CMB anisotropy is often thought to have been generated by quantum mechanical fluctuations in a weakly 
coupled scalar field during an early epoch of inflation and thus would be a realization of a spatially stationary Gaussian 
random process (see, e.g., Ratra 1985; Fischler, Ratra, & Susskind 1985). Measurements lend fairly strong support 
to this Gaussianity assumption (see, e.g., Park et al. 2001; Shandarin et al. 2002; Santos et al. 2002; Polenta et 
al. 2002). See Park, Park, & Ratra (2002) for the effects of foreground contamination on Gaussianity tests based on 
anticipated MAP data. 
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of possibilities. These templates were also used in the preliminary foreground analysis of Coble et 
al. (1999). 

Each of the Galactic emission maps have been converted into a template for cross correlation 
with the PyV data by simulating the PyV observing strategy on it, taking account of the asymmetric 
beam. Since a chopper synchronous offset and a ground shield offset were removed from the data, we 
account for this by adding the chopper synchronous offset and ground shield constraint matrices to 
the noise matrix and marginalizing over them in the analysis. The CMB theory covariance matrix is 
modelled using a spatially-flat cosmological constant dominated CDM model with non-relativistic 
matter density parameter CIq = 0.3, scaled baryonic matter density parameter S1b/i 2 = 0.021 (here 
h is the present value of the Hubble parameter in units of 100 km s _1 Mpc -1 ), and age to = 14 
Gyr (Ratra et al. 1999a) with the CMB anisotropy normalized to a quadrupole moment amplitude 
Qrms-PS = 20 iiK. 



2. Correlations 
2.1. Method 

We follow the general method outlined for example in Hamilton &: Ganga (2001) and Mukherjee 
et al. (2002a), doing a multi-modulation analysis here rather than the multi- (frequency) channel 
analysis discussed in these papers. The method assumes that the data are a linear combination of 
CMB anisotropies and foreground components, 

y = aX + xcmb + n. (1) 

Here y is a n mo d x N element vector containing the data, n is the corresponding noise vector, 
XCMB is the CMB signal, and X is a (n mo< j x N) x n mo d element matrix containing the simulated 
foreground template, N being the number of data points per modulation and n mo( j being the number 
of modulations. In a given column X contains mostly zeros except for the rows corresponding to 
that modulation, where it contains the simulated template signal. The vector a contains n mo( i 
elements that represent the amplitude of the correlated foreground signal. If the noise and CMB 
anisotropies are uncorrelated Gaussian distributed variables with zero mean, minimizing \ 2 leads 
to a best fit estimate of 

a = [X T C- 1 X]- 1 X T C' 1 y, (2) 

where C is the (n mo d x N) x (n mo d x N) element total covariance matrix (sum of the theory covariance 
matrix that models the CMB signal, the noise covariance matrix, and any constraint matrices). The 
vector d then contains the best estimate for the correlation slopes for the corresponding template 
for all modulations, given all information about inter-modulation correlations. The matrix 



^ = {a 2 )-{a) 2 = [X T C- l X}-\ 



(3) 
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is its covariance matrix. The rms amplitudes of temperature fluctuations in the data that re- 
sults from the correlation is AT = (a ± 5a)af ore , where <7f orc contains the rms deviations of the 
corresponding foreground template in the different modulations. 

The method assumes that our foreground emission maps are good enough to model the fore- 
grounds in the data accurately in all the different modulations, at a frequency different from that of 
the original foreground emission map, and that eq. (1) explains all the structure in all modulations 
of the data. If we further assume that the ratio of the signal in the data and that in the foreground 
template is the same for all the modulations, then a net correlation slope can be found using 

Tbtal(£ _1 a) 

° ~ Total(S-i) • 1 ' 

Here Total denotes the sum of all elements of a matrix or vector. This is thus a weighted average 
taking account of correlations between the a values of the different modulations. We may choose not 
to treat all modulations the same (see discussion in § 2.3). The error bar on a is y 7 (Tota^E -1 )) -1 . 



2.2. Results 

The correlation slopes obtained from a complete inter-modulation analysis, using all the PyV 
data points, for each foreground emission template individually, are given in Table 1 (top panel). 
We do not find any significant correlation. The net correlation slopes (i.e., the weighted mean of 
the correlation slopes given in the table, taking cross-modulation correlations into account) are not 
significant, with two standard deviation upper limits of 34 /iK(MJy/sr) -1 , 110 /iK(MJy/sr) _1 , 156 
/jK/fiK, and 897 /zK(MJy/sr) _1 , for the 100 /xm, the 12 (im, the PMN, and the PMNO templates, 
respectively. These limits on 100 fxm. and PMNO correlation slopes compare well with those of Table 
2 of Coble et al. (1999). (There are some errors in the PMN correlation slopes of Table 2 and the 
rms values in Table 3 of Coble et al. 1999.) The upper limits obtained from using all modulations 
of the data together are given in Table 3. 

The results of repeating this analysis for just the SP94 overlap region are also given in Table 
1 (bottom panel). 10 The uncertainties in the estimated correlation slopes are higher here, with 
the number of data points down by a factor of almost 8 (the time spent observing this region was 
less than 1:8 of the total observing time), but the level of associated temperature fluctuations are 
relatively higher, at least for the 100 ^m template. As seen from Tables 2 and 3, the rms of the 
data and the 100 fim template are higher in this patch of sky. Again, no significant correlations 
are found. The weighted correlation coefficient between foreground emission templates is given by 
Er. 1 (E^ 1 S^ 1 )-°- 5 (de Oliveira-Costa et al. 1999). The 100 and 12 fim templates and the PMN 



1 'An outlier that affects 6 data points per feed per modulation is seen in the 12 /im template, hence these data 
points are ignored whenever this region is being analysed. Outliers are harder to pick out over larger regions of sky 
but they also affect the result less, so nothing is removed when analysing the full data set. 
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and PMNO templates are found to have significant weighted correlation coefficients in this region 
of sky, with the correlation reducing somewhat, from 0.7 to 0.5, with increasing modulation for 
the two dust templates, while it remains steady at about 0.85 for the two point source templates. 
(These templates are much less like each other when all the PyV fields are considered.) Yet when 
all modulations of the data are analysed together, joint template fits do not detect any correlations. 
The net correlation slopes obtained from a complete inter-modulation analysis of the SP94 overlap 
region are not significant, with 2 a upper limits of 69 /iK(MJy/sr) _1 , 1027 /iK(MJy/sr) -1 , 603 
fiK/fiK, and 2198 //K(MJy/sr) -1 , for the 100 fim, the 12 /mi, the PMN, and the PMNO templates, 
respectively. The upper limits obtained from using all modulations of the data together are given 
in Table 3. 

2.3. Low I Modulations 

As seen from Table 2, the signal in the foreground templates falls more steeply with increasing 
modulation number than does the data. So the case here, of fitting for several modulations simulta- 
neously, is somewhat different from the case of fitting multi-frequency data: if a certain foreground 
is present in the data, picking it out in the higher modulations will be harder. But the cumulative 
effect of the other modulations is significant on the estimated correlation slope for any modulation; 
Figure 2 of Coble et al. (1999) shows how much the CMB signal in different modulations overlap. 
Hence simultaneously fitting for the correlation slopes in all the modulations of the data may not 
be the most appropriate thing to do. It is also important to note that the foregrounds may not 
have been accurately modelled in the higher modulations because of insufficient resolution (given 
the dense sampling in the data). There may also be different kinds of unmodelled or incorrectly 
modelled foregrounds/errors in different modulations. Hence it is useful to also look at the results 
from analysing say the first three modulations together, because across these modulations the data 
and the foreground signal seem to be roughly similar as regards rms values, or each modulation 
individually, at the cost of increased uncertainty in the estimates. 

If we look at only the first modulation in the SP94 overlap region, we find a 1.8 a correlation 
slope of 80 ± 45 ^(MJy/sr)" 1 (36 ± 21 j*K) for the 100 template (Table 4). For the 12 fim 
template the correlation found is just less than a sigma, but here the two dust templates have a 
weighted correlation coefficient of 0.65, and performing a joint correlation of the data with these 
two templates results in a 1.9 a correlation slope of 115 ± 60 /xK(MJy/sr) _1 (52 ± 27 /iK) for the 
100 /im template at the cost of a 1 a negative correlation with the 12 template. This may be 
the signal found in the foregrounds analysis of the SP94 data (Hamilton & Ganga 2001; Mukherjee 
et al. 2002a; for earlier qualitative indications see Ganga et al. 1997; Ratra et al. 1999b). While 
there a detection was aided by the presence of 7 frequency channels, here the error bars for the 
individual modulation correlations are large. This could be just a chance correlation on the other 
hand, and the probability of that is reflected in the significance of the result. It might also be 
relevant to note here that the SP94 data was only one dimensional while here the PyV 'sa', 'sb', 
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and 'sc' data are essentially two dimensional. 

If we look at only the second modulation by itself, the correlation slopes are as shown in Table 
4, and even though the correlation coefficient between the 12 and 100 fim templates is 0.70 nothing 
is gained by a joint fit this time. 

In the first two modulations together some correlation with the 100 /xm template is detected 
consistent with the above estimates. When the first three modulations are correlated simultane- 
ously, a 1.6 a correlation slope of 58 ± 36 /xK(MJy/sr) _1 (25 ± 16 fiK) is found with the 100 /jm 
template in the first modulation and a 1 a 710 ± 695 / t/K(MJy/sr) _1 (14 ± 14 /jK.) correlation is 
found with the 12 /im template in the second modulation. The correlation coefficient between the 
templates is high and correlations of similar significance are found in the same modulations when 
the two templates are analysed jointly. 

Regarding the point source templates, nothing significant shows up in an individual template 
analysis (Table 4). When analysing the first two modulations together the two point source tem- 
plates have a correlation coefficient of 0.8 and a correlation slope of 6555 ± 4994 ^K(MJy/sr) _1 
(28 ± 21 fiK) at 1.3 a is found with the PMN0 template in modulation 1 at the cost of a 1 a 
negative correlation with the PMN template in the same modulation. When analysing the first 
three modulations together the point source templates have a high correlation coefficient and a 1.2 
a correlation of 5625 ± 4498 /xK(MJy/sr) _1 (24 ± 19 fiK) is found with the PMN0 template in 
modulation 1 at the cost of a —1 a correlation with the PMN template in the same modulation. 

When all the PyV fields are taken together (690 in each modulation), a 1 a correlation with 
the 100 ^m template shows up in the first modulation (when the first modulation is analysed by 
itself, or jointly with the second, or jointly with the second and third), and aid correlation with 
the PMN0 template shows up in the first modulation when it is analysed together with the second 
modulation and jointly with the PMN template. The uncertainty in the correlation slopes with 
the point source templates is lowest in the third and fourth modulations, but nothing significant 
shows up when these modulations are analysed separately. Hence consistent correlations (in the 
low-Z modulations and at low significance) are seen even when both the PyV regions are analysed 
together. 

2.4. Summary 

When all the modulations of the data are analysed together, no significant correlations are 
found, in the whole PyV data set or in the SP94 overlap region. These results are summarized in 
Table 3 and plotted in Figure 1. 

We have discussed the motivations for looking at individual modulations separately, and we 
find some indication of foregrounds in the low / modulations. Correlation slopes for some of these 
low I modulations are summarized in Table 4. We see that greater than 1 a correlations show up in 
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the SP94 overlap region with the 100 /tm template in particular, but also with the 12 /zm template 
and the PMNO template in some modulations in a joint modulation analysis of a few modulations. 
And consistent correlations are seen when all the PyV data fields are analysed together. 

We note that according to the foreground model discussed for example in Mukherjee et al. 
(2002a), if the 100 //m correlations are from both free- free emission and spinning dust emission, 
then by 40 GHz the contribution from spinning dust emission (as traced by the 12 /im template) 
may have again dropped below that from free-free (the frequency at which spinning dust emission 
peaks depends on the details of the spinning dust emission model and this can also vary from 
region to region), so that we would expect 100 /im correlations to be more significant than 12 /im 
correlations at this frequency (see Figure la). However, given the error bars (the numbers in Table 
4), we do not expect to see significant correlations in the PyV data (see Figure 1). 

3. Conclusion 

Using the method of cross correlating CMB anisotropy data with foreground contaminant 
emission templates, we find little evidence for foreground contamination in the whole PyV region 
when analysing all modulations of the data together. There is however indication of foreground 
contamination in the low-Z modulations of the PyV fields that cover the region scanned earlier 
by the SP94 experiment, consistent with results from the SP94 data. This is a valuable test, 
indicating consistency between results found using data from two different experiments. Given the 
uncertainties, our findings are not inconsistent with the two component dust-correlated (free-free 
and spinning dust) foreground emission model that other data sets tentatively seem to point to 
(Mukherjee et al. 2002a). 

We thank D. Finkbeiner for the 12 /mi dust map and E. Boyce, R. Ekers, E. Ryan- Weber, 
and L. Stavely-Smith for their work on this project. PM, BR, and TS acknowledge support from 
NSF CAREER grant AST-9875031. KC is supported by NSF grant AST-0104465. This work 
was partially carried out at the Infrared Processing and Analysis Center and the Jet Propulsion 
Laboratory of the California Institute of Technology, under a contract with the National Aeronautics 
and Space Administration. 



Table 1. Results of Correlating the PyV Data with Individual Foreground Emission Templates a 



Template 


modi 


mod2 


mod3 


mod4 


mod5 


mod6 


mod7 


mod8 


100 /im 


1 ±20 


-4 ±20 


-20 ± 20 


-6 ±23 


5 ±30 


8 ±52 


-34 ±110 


-172 ±245 




0±4 


-1±3 


-2 ±2 


0±2 


Oil 


Oil 


Oil 


-1 ± 1 


12 /im 


-82 ±131 


-75 ± 123 


-79 ± 120 


-133 ± 126 


-257 ± 140 


-372 i 182 


-228 i 259 


-476 i 516 




-2 ±3 


-2 ±3 


-1±2 


-2 ±2 


-3± 1 


-3i 1 


-1 i 1 


-1 i 1 


PMN 


137 ± 123 


66 ± 108 


-14 ± 103 


-68 ± 100 


-158 ± 105 


-173 i 126 


-137 i 177 


8 ± 252 




3±3 


2±3 


0±2 


-1±2 


-2± 1 


-2i 1 


-1 i 1 


Oil 


PMNO 


684 ± 502 


546 ± 408 


237 ± 382 


15 ± 369 


-111 ±385 


-116 i 452 


-44 i 625 


972 i 986 




3 ± 2 


3±2 


1 ±2 


Oil 


Oil 


Oil 


Oil 


1 i 1 


100 /im 


19 ±30 


18 ±30 


5 ±31 


12i35 


30i50 


92i94 


-104 i 266 


792 i 1051 




8± 14 


6 ± 11 


1±7 


2±6 


3±5 


4i4 


-2 ±4 


4i5 


12 pLm. 


116 ±570 


235 ± 559 


-288 ± 561 


-248 ± 647 


-348 ± 868 


-626 i 1672 


-2476 ± 3871 


7207 ± 11467 




2 ± 12 


4± 10 


-4 ±7 


-2 ±6 


-2 ±6 


-2i5 


-3 ±5 


4i6 


PMN 


205 ± 327 


98 ± 291 


70 ± 266 


5 ±276 


78 ± 320 


-123 ±436 


319 ± 801 


862 ± 1717 




5±8 


3±8 


1 ±6 


0±4 


1 ±4 


-1±3 


1 ±4 


3i5 


PMNO 


1871 ± 2030 


947 ± 1722 


636 ± 1607 


267 ± 1645 


409 ± 1958 


-1496 ± 2851 


3495 ± 5538 


9599 ± 11042 




8±8 


5±8 


2±6 


1±5 


1 ±4 


-2 ±4 


3i4 


5i6 



a For each template the correlation slopes (in units of /xK(MJy/sr) _1 for the dust and the PMNO templates and in units of 
/xK//xK for the PMN template) are given in the first row and the corresponding temperature fluctuations (in /xK) are given in the 
second row. In the top panel we have used all the data points, and in the bottom panel we use only the fields in the SP94 overlap 
region. Cross-modulation correlations are accounted for and the 8 modulations represent decreasing angular scales covering / ~ 40 
to ~ 260. 
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Table 2. Rms Values of the Signal in the PyV Data and the Templates a 



Modulation 


Data 


100 nm 


12 /j,m 


PMN 


PMN0 




/iK 


MJy/sr 


MJy/sr 


UK 


MJy/sr 


1 


no n 

92. u 


n ooo 
U.zzz 


n noo 
U.Uzz 


n noo 
U.Uzz 


n nn A 
0.004 


z 


n£ i 

95.1 


0.1/0 


n noo 
U.Uzz 


n no^7 
U.Uz r 


n nn£? 
0.00b 


Q 

o 


<y.4 


n i nn 
U.1UU 


u.uiy 


n noo 
U.Uzz 


n nn^ 
U.UUO 


4 


( l.U 


o.o ro 


n m £ 
U.Ulo 


n m ^7 
U.U1 7 


n nn /i 
0.004 


r 




f o.o 


U.U4 f 


n m n 
U.U1U 


n ni q 
U.Ulo 


n nnQ 
U.UUo 


6 


77.4 


0.023 


0.007 


0.010 


0.002 


7 


69.8 


0.011 


0.006 


0.007 


0.002 


8 


65.0 


0.006 


0.003 


0.005 


0.001 


1 


123.0 


0.455 


0.022 


0.025 


0.004 


2 


129.9 


0.356 


0.018 


0.028 


0.005 


3 


104.9 


0.221 


0.013 


0.022 


0.004 


4 


111.1 


0.180 


0.010 


0.016 


0.003 


5 


113.6 


0.103 


0.007 


0.012 


0.002 


6 


111.5 


0.042 


0.003 


0.008 


0.001 


7 


95.7 


0.015 


0.001 


0.004 


0.001 


8 


89.33 


0.005 


0.001 


0.003 


0.0005 



a The top panel is for the whole PyV region with 690 fields 
per modulation. The bottom panel is for the SP94 overlap 
region with 90 fields per modulation. 



Table 3. Correlation Slopes 



Fields 


Modulations 


100 fim 


12 fim 


PMN 


PMN0 






^(MJy/sr)- 1 


^(MJy/sr)- 1 




^(MJy/sr)" 1 


all 


all 


-4± 19 


-120 ± 115 


-30 ± 93 


209 ± 344 


SP94 overlap 


all 


13 ±28 


-21 ±524 


85 ± 259 


904 ± 1551 
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Table 4. Correlation Slopes for Low I Modulations 3, 



Fields 


Modulations 


100 fim 


12 /mi 


PMN 


PMN0 






^(MJy/sr)- 1 


^(MJy/sr)- 1 




^(MJy/sr)- 1 


all 


1st 3 mods 


7 ±22 


-155 ± 141 


-59 ± 131 


203 ± 553 


all 


1st 2 mods 


16 ±24 


-70 ± 162 


-77 ± 154 


416 ± 658 


all 


1st mod 


29 ±29 


97 ± 223 


-161 ±213 


213 ± 903 


all 


2nd mod 


29 ±33 


-27 ±202 


-53 ± 190 


271 ± 779 


SP94 overlap 


1st 3 mods 


33 ±32 


388 ± 590 


-129 ±332 


-283 ± 1924 


SP94 overlap 


1st 2 mods 


50 ±36 


522 ± 674 


-6 ± 387 


736 ± 2215 


SP94 overlap 


1st mod 


80 ±45 


406 ± 875 


190 ± 563 


2131 ± 3401 


SP94 overlap 


2nd mod 


49 ±47 


854 ± 889 


65 ± 473 


-158 ±2805 



a Boldface type indicates 1 a detection. 
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10 20 50 100 

Frequency (GHz) 

Fig. 1. — a) Correlation slopes between CMB anisotropy data and the dust 100 template. Two 
standard deviation upper limits obtained for PyV from analysing all the fields and from analysing 
the SP94 overlap region only are both shown at ~ 40 GHz. b) Correlations with the 12 /xm template. 
Results are also shown for other experiments: Tenerife (open stars, de Oliveira-Costa et al. 2001; 
filled stars, Mukherjee et al. 2001), OVRO (filled triangles), 19GHz (open triangles), SP94 (open 
circles), Saskatoon (open square), Boomerang (filled rectangle), and DMR (filled circles). The 
solid curves are representative of a spinning dust spectrum, the dashed and the dot-dashed lines 
represent free-free and vibrational dust emission spectra, respectively. 



